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Abstract. We review the cosmological implications of a relic population of pseudo Nambu- 
Goldstone bosons (pNGB) with an anomalous coupling to two photons, often called axion-like 
particles (ALPs). We establish constraints on the pNGB mass and two-photon coupling by 
considering big bang nucleosynthesis, the physics of the cosmic microwave background, and 
the diffuse photon background. The bounds from WMAP7 and other large-scale-structure 
data on the effective number of neutrino species can be stronger than the traditional bounds 
from the primordial helium abundance. These bounds, together with those from primordial 
deuterium abundance, constitute the most stringent probes of early decays. 
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1 Overview 

The spontaneous symmetry breaking (SSB) of exact global continuous symmetries produces 
massless Nambu-Goldstone bosons. If the symmetries are only approximate, the puta- 
tively massless bosons acquire small masses and are called pseudo Nambu-Goldstone bosons 
(pNGB). The axion is a well studied example of the latter: it arises from the SSB of the 
Peccei-Quinn U(1)pq axial symmetry, postulated to solve the strong CP problem [1, 2]. This 
symmetry is explicitly broken only by the colour anomaly. The axion potential becomes 
~ Aq CD cos(a// n ), where a is the axion field, f a is the PQ SSB scale up to an integer and 
Aqcd is the confinement scale of quantum chromodynamics. As a consequence the axion 
picks up a mass 1 m a ~ Aq CD // a . At energies much below f a , the axion field always enters 
into the effective Lagrangian in the combination a/f a so its interactions with standard model 
(SM) particles are always suppressed by the high energy scale f a . 

The axion can be generalised to a generic pNGB, which in this paper we call an axion- 
like particle (ALP) and denote by <j>. If the dynamics explicitly breaking the associated global 
continuous symmetry have a characteristic scale A, the ALP mass is rriA ~ A 2 //^,, where 
is the scale at which the ALP symmetry is spontaneously broken. Interactions of ALPs with 
SM particles are also suppressed by this scale f^. The phenomenology of the SM requires A 
to be related to physics beyond the electroweak scale, i.e. A > TeV (which implies S> m a 
for = f a ), or to belong to a hidden sector. Axions, ALPs and other low mass particles 
postulated as a bridge between the hidden sector and the SM are typical examples of weakly 

In nature there is another pNGB, the rj' meson (related to a common axial phase shift transformation of 
the light quarks) with which the axion mixes and through it, with the rest of pseudoscalar mesons. At the 
end of the day, the mass of the axion further reduces to m a ~ m^/^-// a . 



- 1 - 



interacting slim particles (WISPs), which are subjects of growing interest in the low-energy 
high-precision frontier of particle physics [3]. 

Axions and ALPs are subject to very strong constraints from stellar evolution and 
cosmology. The more stringent astrophysical bounds rely on the production of these exotic 
particles in the hot and dense stellar interiors and their subsequent escape contributing 
directly to stellar energy losses and therefore to the consumption of nuclear fuel. Normally 
the weaker the ALP interacts, the smaller the ALP production and the impact on stars. 

In cosmology, the situation is reversed. Assuming sufficiently high initial temperature 
of the universe, a thermal population of ALPs is created and decouples from the SM thermal 
bath at some point. The impact on the cosmological observables at our disposal — big bang 
nucleosynthesis (BBN), anisotropies in the cosmic microwave background (CMB), etc. — 
mainly depends on when this relic population decays (if it decays at all) . The decay lifetime 
is proportional to f^/rr&, so the weaker the coupling, the slower the decay, and the greater 
the relevance to late cosmology. Considering this, cosmological arguments are complementary 
to stellar evolution constraints and laboratory experiments. 

In a recent paper [4], we studied the constraints that current cosmology places on 
cosmologically unstable axions (those with masses in the range of a few eV to MeV). In this 
paper we extend this study to a general ALP. In [4] we found that the phenomenology is 
basically determined by the anomaly-driven two-photon coupling, which for an ALP we write 
as 

C Hl = -^F^F^<P (1.1) 

where F^ v and F^ v are the electromagnetic field tensor and its dual. The coupling strength 
g ~ a/(2irf^) has inverse mass dimension and provides the ALP with an efficient thermal- 
ization mechanism in the early universe, i.e. the Primakoff process 7 + <? — > (ft + q (q stands 
for any charged particle), and a prominent decay channel, i.e. (ft — > 77. 

For most of the paper, we focus on this coupling, which we take as the defining char- 
acteristic of an ALP, together with its mass. The effects produced by other couplings are 
finally discussed in Sec. 5. 

Alongside the thermal ALP population, a non-thermal population of ALPs can be cre- 
ated by the realignment mechanism just as in the axion case [5]. The consequences of this 
population will be presented in a separate publication [6] . 

A number of excellent papers have considered the cosmological implications of decaying 
massive relics, see for instance [7-9]. Our work focuses on ALPs and deals with particles 
in an intermediate mass range (keV-GeV), which is often not considered. As an additional 
complication, in a part of parameter space the inverse decay 77 — > (ft is also effective. ALPs 
might not completely decouple nor be in thermal equilibrium so we must follow the evolution 
of the ALP phase space distribution. This paper updates and complements the previous 
work of Masso and Toldra on a spinless particle coupled to photons [10, 11]. 

Our conclusions are presented in Figs. 1 and 2, which shows the constraints on ALPs 

—3/2 

coupled to photons. Cosmological bounds in Fig. 1 have a characteristic slope g oc 
because they critically depend on the ALP lifetime 

(1.2) 

and have a milder dependence on other parameters. Ordered by decreasing lifetimes, the 
excluded regions are: 
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• DM — if ALPs are cosmologically stable and behave as dark matter (DM) they should 
not exceed the DM fraction measured by WMAP; 

• Optical, X-Rays, 7-Rays — photons produced in ALP decays inside galaxies would 
show up as a peak in galactic spectra that must not exceed the known backgrounds; 

• EBL — photons produced in ALP decays when the universe is transparent must not 
exceed the extragalactic background light (EBL); 

• Xi on — the ionization of primordial hydrogen caused by the decay photons must not 
contribute significantly to the optical depth after recombination; 

• CMB y, /x — if the decay happens when the universe is opaque, the decay photons must 
not cause spectral distortions in the CMB spectrum that cannot be fully rethermalized; 

• EM, Hadr showers — the decay of high mass ALPs produces electromagnetic and 
hadronic showers that must not spoil the agreement of big bang nucleosynthesis with 
observations of primordial nuclei; 

• 4 He, D — the ALP decays produce photons (entropy) that dilute the baryon and 
neutrino densities, whose values affect the outcome of BBN, in particular the deuterium 
and 4 He yields. Again, this dilution should not compromise BBN; 

• N e ff — the neutrino density must not disagree with the value measured by WMAP 
and other large-scale-structure probes. Currently, data points to a number of effective 
neutrinos iV e g greater than 3, which is disfavoured in the decaying ALP cosmology. 

These bounds are complemented at low masses and large couplings by stellar-evolution 
arguments and laboratory searches. The most relevant astrophysical bounds come from star 
counts in the colour-magnitude diagrams of globular clusters, in particular through the effect 
on the evolution of horizontal branch stars. This constraint overlaps significantly with the 
cosmological bounds, so we have computed precisely its turn off at large ALP masses (region 
labelled HB in Figs. 1 and 2). The duration of the neutrino pulse from supernova 1987a 
can also constrain ALPs (SN) [10]. However this bound relies on very sparse data and 
insufficient understanding of supernova dynamics and ALP emission from a nuclear-density 
environment, which makes it somehow unreliable. This is a pity because the temperature 
and density of the SN core at the time of the collapse are among the highest we can find 
in stellar environments, so the Boltzmann suppression of the bounds appears only at larger 
masses [12]. 

Much stronger couplings are tested by direct searches in laboratories (for a review 
see [3, 13]), however the sky is the fundamental tool to test weakly coupled ALPs. 

We have divided this paper in three main blocks. In Sec. 2 we discuss the relic abun- 
dance of ALPs. In Sec. 3 we consider the cosmological implications of ALP decays before 
recombination, in the opaque universe. The bound from the effective number of neutrinos 
is the most relevant, since it covers most of the parameter space and constitutes the most 
relevant original contribution of this paper to the ALP parameter space. We also describe 
the implications on BBN. 

We describe the signatures and bounds on ALPs decaying after recombination in Sec. 4. 
We present bounds from the non observation of the decay photons and the ionization history 
of the universe and the possible signatures of ALP decay in the CMB spectrum. Finally, in 
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Sec. 5 we briefly analyse how our bounds change with the addition other couplings besides 
the two-photon one. The computation of the high mass boundary of the HB constraint is 
discussed in the Appendix. The summary plots are Figs. 1 and 2 and our conclusions are 
drawn in Sec. 6. 

2 Establishment and decay of a thermal ALP population 

The two-photon coupling, Eq. 1.1, allows the establishment of a thermal population of ALPs 
in the early universe through the Primakoff process. The rate of ALP production due to 
scattering on relativistic electrons was computed in [14] to be 

r '=if r3 H^) + °H' <"> 

where m 7 = eT/3 is the plasmon mass in an electron-positron plasma and T the temperature. 
Taking the rate to be proportional to the number density of electrons, n e = 3£(3)T 3 /it 2 , we 
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Figure 2. ALP bounds in the — r parameter space. The labeling is described in Sec. 1. 



can generalise to a multicomponent plasma, 



gg 2 k 2 . 
r„ ~ i log 



36C(3) 



T 2 

2 



in 



+ 0.82 n 



(2.2) 



where n q is the effective number density of charged particles, n q = ^ Qf ni = (C(3) /7r 2 )g q (T)T 3 , 
where Qi is the charge of i-th particle species, and the parameter g q (T) represents the ef- 
fective number of relativistic charged degrees of freedom. Also the plasmon mass has to be 

1/2 

corrected by a factor m 7 oc g q when more charged species are present. 

The Primakoff process becomes inefficient when T q becomes smaller than the Hubble 
expansion parameter. In a radiation dominated universe, this is H = 1.66^/5* (T)T 2 /m P i 
where mpi is the Planck mass and g*{T) is the effective number of relativistic degrees of 
freedom, which like g q {T) is a smooth function of T. We follow the definitions and notations 
of [15] for all the usual cosmological quantities. The T dependences of g q and g* are taken 
from the recent analysis of [16]. The ALP bath decouples at a temperature very sensitive to 

Tfo B _JL_VE _ 123 v^ ( m-^-'V GeV. (2.3) 



atg "mpi g q 



9g 
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For values of g < 2 x 10 -9 GeV, ALP interactions freeze out at temperatures above the 
electroweak scale (Eew ~ 250 GeV), where the particle content of the plasma is relatively 
speculative. For instance in the MSSM scenario, above the SUSY breaking energy scale we 
have <?* = 228.75, while the SM alone provides only 106.75 relativistic degrees of freedom. 
We warn the reader that for g < 10" 17 GeV we require a most likely meaningless freeze out 
temperature above the Planck scale (mpi = 1.22 x 10 19 GeV). We drew two dotted black 
lines for Tf Q = Eew an d for Tf Q = mpi in Figs. 1 and 2. 

After the freeze out, the number density of ALPs per comoving volume is conserved. 
Also its ratio with the comoving entropy density = N^/S is conserved. Therefore at any 
later time (before the ALP decay), the ALP density will be given by 

n^(T)=Y 4> (Tf )s(T). (2.4) 

Given the standard particle content that we are assuming, there is a minimum ALP yield 
given by the value of g*(T{ > EW scale) 

^>^^- 0.005 9 ,(T). (2.5) 
n 7 2 10b. 75 

Cosmologically stable ALPs must not exceed the measured abundance of DM, therefore 

n h 2 = P± h 2 = HHJH h 2 < nBuh 2 = Q U (2 6) 

Pc Pc 

where p c = 3H 2 mp^/8ir is the critical density and h 2 is the normalised Hubble expansion 
rate. ALPs with a mass = 154 eV would account for all the dark matter of the universe 
but larger masses are excluded. This exclusion is depicted in orange in our summary plots 
Figs. 1 and 2 and labelled as DM. The discovery of new degrees of freedom (dof) above 
the EW scale would relax this bound, which is linearly sensitive to <?*(7f ), by a factor 
(106.75 + new dof)/106.75. This means that 0(100) of them are needed for a sizeable 
change. 

Let us now consider the dynamics of the ALP decay. The guiding quantity is the ratio of 
the decay rate over the Hubble parameter evaluated at 3T ~ (when ALPs would become 
nonrelativistic if in thermal equilibrium) disregarding the ALP contribution 

r 7 ^ 3 - 3 m * ( a y (27) 

H(m^/3) Vff* ("V 3 ) MeV V lO" 7 GeV" 1 / V ' 7 

If this ratio is large, the condition T^/H ~ 1 is satisfied at a temperature higher than 
since H grows with T. In these conditions, the inverse decay process, 77—7-0 (photon 
coalescence) is effective and the ALP population regains thermal equilibrium [4]. Note that 
during this rethermalization process, a small amount of entropy will be generated. The 
population of ALPs will be kept in thermal equilibrium with photons by decays and inverse 
decays. When the temperature decreases to T ~ m^/3, the thermal abundance starts to 
become Boltzmann suppressed and axions disappear from the bath in thermal equilibrium, 
very much like electrons and positrons annihilate at T ~ m e /3. 

On the other hand, if the ratio is small, T~/H ~ 1 will happen when T <C and ALPs 
are very nonrelativistic. In this case, inverse decays are Boltzmann suppressed and the ALP 
decay is not perturbed by the thermal bath of standard particles. It is also possible that 
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the energy density in ALPs, m^n^, becomes larger than that of standard model particles, 
(7r 2 /30)<7*(T)T 4 , and the ALP decays create a considerable amount of entropy. This happens 
if the decay temperature T d satisfies 

m 4> > 71-4 9*(Tio) _ nr7 9*(T{ ) 



T d ~ 30C(3) g*(T d ) ' g*(T d 

If ALPs do not dominate the universe energy budget when they decay T d can be esti- 
mated by T 7 ~ H(T d ), which gives 

T d - 1/ °; 6 (—^ T ) f^l^MeV. (2.9) 

gJ (T d ) V10- 7 GeV-V VMeV J V ; 

If the ALP energy density does dominate, the decay temperature is instead 

Td °A ( i V /3 f^V /3 MeV (2 10) 

which is typically larger than the previous case (a matter dominated universe expands more 
slowly than a radiation dominated one). If the universe becomes radiation dominated after 
the decay, the temperature in Eq. 2.9 gives the correct order of magnitude for the reheating 
temperature. We evaluate the observable consequences of early ALP decay in the next 
sections. 



3 Early ALP decays 
3.1 Neutrino dilution 

The ALP population decay produces photons whose energy and entropy tends to be reshuf- 
fled between the different species that are present and active in the thermal bath at that 
moment. The relevant temperatures in the context of this paper are around and below the 
neutrino decoupling temperature, T ~ a few MeV. The active species to consider are then 
photons, electrons and neutrinos. Photons thermalize very fast with electrons, the particles 
most tightly coupled to neutrinos. Weak reactions between electrons and neutrinos heat 
the neutrino bath to keep track with the electron/photon temperature changes. The most 
relevant of these energy redistribution processes is e + e~ — > vv. Since scattering processes 
such as e^v — > e v cannot change the neutrino number and are less effective in transferring 
energy to the neutrino bath, we neglect them 2 . Therefore, if ALP decay happens after the 
freeze-out of the e + e~ f-> vv reactions, the electromagnetic energy and the entropy would 
not be shared with the neutrino bath, which then appears to have less energy than in the 
standard case. The energy flow from ALPs to neutrinos can then be modelled by a set of 
Boltzmann equations for comoving energies, defined as 

Xi = Pi R\ (3.1) 

2 Assume that electrons have a larger temperature. Then the speed of energy transfer per unit volume 
in scattering processes will be oc {5E)T^T* with (5E) a thermal-averaged energy transfer per scattering 
(~ T e — T„); while for annihilations it will be ~ T e (T^ — ), which is much more sensitive to nonequilibrium 
situations. 
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where R is the cosmic scale factor 3 , and the ALP phase space distribution function f{k&) as 
a function of the comoving momentum, kj, = Rh^, 



-f(k rl> ) = -(c J + c q )(f-D, 

j- t {X, + X e ) = 3H5 Pe R* + |^^(C 7 + C q )(f - / cq ) 

+ T ^{c e (xl-xzf)+c» T (xl T -x Z 
Tt x„ e -- w c e {x Ue -x Ue ) , 



d -R=iJ J^r- [X 1 +X e + X Ue + X^ T + p^) , 



H" Ry 3rripj 



(3.2a) 

(3.2b) 
(3.2c) 
(3.2d) 

(3.2e) 



where = y &| + m 2 is the ALP energy and Cj^ = uj^R. The collision terms for the decay 
and Primakoff processes are [4] 



( 



nr 



4m 7 m 4> 



lb 



1 + 



2T 1 
1 + — log - 
1 



e -(w -^)/2T 
2 



1 » 



(4^(m e + 3r)) 
m 2 (m 2 + (m e + 3T) 2 ) 



(3.3) 
(3.4) 



The equations 3.2 describe the evolution of the comoving energy density stored in 7 
together with e^, the three species of i/s, the ALPs and the cosmic scale factor R. The 
electron-neutrino energy exchange rate depends on the neutrino flavour because of the ab- 
sence of charge current interactions for p and r flavours. At the low temperatures of interest 
F eu = GpT^ with Gf the Fermi constant, and C e ~ 0.68 and Cn T ~ 0.15 which follows from 
the appropriated thermally averaged cross section. The factor 5p e = p e — p e /3 (p e and p e are 
the pressure and energy density of e ± ) accounts for the comoving energy density gain as elec- 
trons become increasingly nonrelativistic. We assume that neutrinos always have a thermal 
distribution, determined only by an effective temperature, which should be a reasonable first 
approximation and accurate enough for our purposes. We also neglect the energy reshuffling 
between different neutrino species, which does not influence the total neutrino density at 
leading order. The initial conditions are specified at T > MeV by having all species at a 
common temperature and the ALP number density given by Eq. 2.4. 

For values g < 10 _7 GeV _1 , the Primakoff process is decoupled in the temperature range 
of interest and can be neglected. If the inverse decay is also negligible, we can integrate the 
ALP phase space distribution explicitly and, instead of using Eq. 3.2a, directly compute the 
evolution of the number density 



d 



— (n^R 3 ) = — r 7 n^i? , 



(3.5) 



R has dimension of [energy] 1 and so the X's are dimensionless. 
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Figure 3. Contour plots of the number of effective neutrinos N c g as a function of the ALP mass 
and lifetime (left) and of the coupling parameter g and lifetime r (right). ALP cosmologies leading 
to JVoff < 2.11 can be safely excluded. 



recovering the exponential decay law oc e 7 *. The integral in Eq. 3.2b is then 

7^4^(C 7 + C q )(f - f*) ~ m^n^R 4 . (3.6) 

The final neutrino energy density is usually quoted as the effective number of standard 
neutrinos 

^ = ^W' <3 ' 7) 

8Uli ^7 

We have scanned the ALP parameter space and present our results for jV e ff in the g — r 
and m<f) — r planes in Fig. 3. In Fig. 4, we present some illustrative examples depicting the 
evolution of the X's of electrons, neutrinos and ALPs as a function of the temperature. Note 
that in all of them, when ALPs become nonrelativistic, the ratio X^/X^ rises because it 
becomes proportional to m^/T, until the age of the universe becomes comparable with r. 

Fig. 4(a) shows a typical case of very early decay of a massive ALP, when neutrinos are 
still partially coupled to electrons. Here N e g = 2.6, marginally different from the standard 
value of 3. Indeed, even if the ALP energy dominates the universe and the entropy injected 
during the decay is humongous, the reheating temperature is large enough for neutrinos to 
almost fully recover their thermal abundance. If the ALPs decay earlier the outcome is 
identical to standard cosmology since neutrinos will regain completely their thermal abun- 
dance. In this and similar cases, the final value of N e g is related not to the total entropy, 
but only to that part injected after the freeze out of the neutrino-electron interactions. The 
neutrino dilution is sensitive mainly to the ALP lifetime and not to or g individually. 
This is precisely what we find in the lower right corner of Fig. 3(a), and correspondingly to 
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the lower left of Fig. 3(b). Around T ~ m e , electrons and positrons become nonrelativistic 
and annihilate, heating the photon bath (and baryons) but not the neutrinos, which have 
decoupled. The ratio X u /X 7 therefore decreases. In this period, the photon temperature 
increases with respect to the neutrino temperature by the standard factor (4/1 l) 1 / 3 due to 
entropy conservation 

(T 7 /r^)' = (4/11)^(^/1],). (3.8) 

The picture changes if we consider later decays. At a late enough point, the decays 
proceed when the neutrinos have already decoupled; photons and electrons get all the ALP 
energy, as shown in Fig. 4(b). The neutrino dilution is computable in this case. The ratio of 
the final and initial comoving entropies of the photon/electron bath is [15] 



^ = 1.83 (ff 1/3 ) 3/4 mY(t> 
Si * s ^mpif. 



^ = 1.83^)""-^. (3.9) 



1/3 

where (gj s ), the time average of g*s during the decay, is an 0(1) factor. The temperature 
of the electromagnetic bath increases with respect to the neutrino one by a factor (Sy/S'i) 1 / 3 
making A e g- = 3(S , j/5j) _4//3 because also the electron entropy ends up in photons. The 
neutrino energy density is therefore strongly diluted by the energy gain of photons plus 
electrons. Note that this is mainly a function of m^j '-^/Fy oc (g^/m^) -1 , which produces the 
characteristic slope of the isocontours at long ALP lifetimes r oc m~ 2 for Fig. 3(a) and r oc g 4 
for Fig. 3(b). For g < 10 -9 GeV -1 , this is the only dependence on g, since is constant. 
Another example of late ALP decay, but with a smaller mass, is shown in Fig. 4(c). In this 
case we observe first the annihilation, which heats photons with respect to neutrinos. A 
sizeable neutrino dilution is observable after the ALP decay. 

The relic abundance of ALPs grows for bigger values of g. In particular, around g ~ 
10~ 7 5 GeV~ 1 it suffers an abrupt increase due to the sizeable decrease of g* during the QCD 
confining phase transition. The isocontours in Fig. 3 have sharp features in this g range 
(in Fig. 3(a) the g dependence is hidden in r). ALPs decoupling at smaller temperatures, 
i.e. with bigger g values, do not get their abundance diluted by the QCD degrees of freedom: 
being more abundant, they produce more entropy when they decay. 

Finally, in Fig. 4(d) we show an example for which inverse decays are relevant. As 
the temperature drops, we observe a first decrease of the ALP energy density due to the 
electrons heating the photon bath. The inverse decay channel opens around T ~ 70 keV and 
helps ALPs to regain equilibrium before disappearing at T ~ m^. During rethermalization, 
the photon energy decreases, which can be seen as a slight rise in X u /X^ (neutrinos are 
decoupled by then). Neglecting the entropy gain in the mixture process, which is typically 
small, entropy conservation gives N c q = 3(11/13) 4 / 3 ~ 2.4 [4]. Due to this mechanism, in the 
small-mass and short-lifetime region of the parameter space the prediction is N e g ~ 2.4. If 
m <t> ~ a f ew MeV, the disappearance from the thermal bath happens when neutrinos are still 
coupled, so iV e flf approaches 3. Since the disappearance of ALPs in local thermal equilibrium 
is governed only by the ALP mass, the isocontours of N e ^ exactly follow the isocontours of 

From the observational point of view, the presence of neutrinos modifies structure for- 
mation and leaves a trace in the anisotropies of the CMB [15, 18]. The two main effects are 
a shift in the redshift of matter-radiation equality and the decrease of the density contrast 
at small scales caused by the neutrino free-streaming [19-22]. In our previous paper, we 
derived constraints from WMAP7 and other cosmological data on N c r assuming a flat prior 
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Figure 4. Examples of the evolution of the comoving energy of all species of neutrinos (solid), ALPs 
(dashed) and electrons (double-dashed), as functions of the temperature of the universe. All energies 
are normalized to one thermal bosonic degree of freedom. 



< N e g < 3. This is the range allowed in our scenario provided that ALPs have decayed 
before matter-radiation equality and the standard analysis of CMB is unperturbed by the 
presence of ALPs. See [4] for details and references. We found 



N cS > { 



2.70 at 68% C.L. 
2.39 at 95% C.L. 
2.11 at 99% C.L. 



(3.10) 



Our numerical calculations exclude the yellow region in Figs. 1 and 2 (labelled iV e fr) at 99% 
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C.L.. Remarkably, the 95% C.L. is just below the value N e g = 2.4, predicted when ALPs 
are in thermal contact through Primakoff and inverse decays. The region disfavoured at the 
95% C.L. is much bigger than at 99% C.L. but is extremely sensitive to the exact value of 
the bound so we prefer to quote the 99% C.L.. As noted in [4], the Planck satellite will not 
improve this bound, if it measures 7V e g = 3 with the predicted error bars. The strength of 
the WMAP7 bound relies on cosmology seeming to prefer N c g > 3 [23]. Of course, if the 
Planck satellite, which is presently taking data, finds convincing evidence for extra radiation, 
much in cosmology will have to be reconsidered besides our ALP limits. 

The mass lower bound of the bound is given by the limit of validity of our assumptions, 
i.e. standard cosmology at temperatures below the standard matter-radiation equality T ~ 1 
eV. From Eq. 2.9 we estimate this to be 



3.2 Big bang nucleosynthesis 

The decays of relic particles can have a considerable effect on the abundances of light elements 
predicted by big bang nucleosynthesis (see [33-35]). The impact depends strongly on the ALP 
mass, in particular whether ALPs are heavy enough to induce electromagnetic or hadronic 
cascades. We discuss these cases separately. 

Small masses 

If the ALP mass is smaller than a few MeV, the decay products cannot induce nuclear 
reactions and their effect on BBN is only indirect. The injected photons (and perhaps a small 
amount of electron/positron pairs) dilute both the neutrino and baryon densities, which in 
turn control the effectiveness of the standard nuclear reactions happening during BBN. The 
impact then depends on whether the decays happen before or after BBN. 

For decays happening after BBN, the injected photons heat the bath, decreasing the 
baryon to photon ratio 4 , r\. Therefore, the value of r\ during BBN is larger than the one 
measured much later, at the CMB epoch by WMAP7 and other large-scale-structure data 
r/CMB = 6 _23 ± 0.17 [23] (i] 10 = 10 10 ??). If BBN proceeds in such a low-photon-density 
environment, the deuterium photo-disintegration reaction is less effective, so that a significant 
amount of deuterium forms earlier, i.e. the deuterium bottleneck opens up earlier. This allows 
the subsequent reactions (burning deuterium into heavier elements) to happen at higher 
baryon densities, where they are more effective. The outcome of this high-r/BBN scenario is 
clear: intermediate nuclei like D or 3 He are more easily consumed, and the final abundance 
of heavier nuclei like Li increases. 

The Helium abundance is determined by the neutron abundance during BBN, since 
essentially all neutrons end up in Helium nuclei. Let us then briefly review the history of 
neutrons. At high temperatures T 3> MeV, electroweak reactions keep protons and neutrons 
in nuclear equilibrium, which predicts n n /n p = e~® vn l T (Q pn ~ m n — m p = 1.293 MeV) for 
the neutron/proton ratio. At T ~ MeV, the reactions freeze out and afterwards neutrons can 

4 The value of r\ can be measured at two different epochs, during BBN (estimated by measuring the 
primordial abundances) and at the CMB epoch (as imprinted in the temperature angular anisotropies). These 
two estimates agree quite well, rtfa™ = 6.23±0.17 [23] and 5.1 < rjf™ < 6.5 at 95% C.L. [2]. This agreement 
can be readily used as a constraint on any entropy injection between the BBN and CMB epochs. However, 
we also want to consider here the cases in which the ALPs decay during BBN, at least partially, where the 
application of this constraint is not sufficiently clear. 
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only decay until the deuterium bottleneck opens up, BBN commences and they get trapped 
in 4 He nuclei. In the late ALP decay scenario, the bottleneck opens earlier and thus less 
neutrons decay, enhancing the final 4 He yield. But this is not the only effect. The presence 
of ALPs makes the universe expand faster, which has two additional implications: a) during 
the freeze-out of the p f-> n conversion reactions it induces an earlier freeze out and therefore 
a larger n abundance and b) it shortens the time between this freeze-out and BBN and 
therefore the amount of neutrons that decay. Thus three mechanisms are responsible for the 
enhancement of the 4 He yield. 

Let us now consider the case where ALPs decay before BBN, i.e. before the opening 
of the deuterium bottleneck. In this case tjbbn = f|CMBi and the main trends mentioned 
before disappear. However, the ALP decays can still modify BBN indirectly, if they happen 
between the freeze-out of weak interactions and BBN. There are three effects that we should 
take into account. First, ALPs are present during the freeze-out of the p f-> n reactions, 
so the n abundance is in principle larger. Second, when ALPs decay they reduce N c g, as 
shown in the previous section. After this event, the universe expands slower than usual and 
therefore: a) more neutrons decay before BBN and b) BBN happens in a low iV e g. 

As in the case of a post-BBN ALP decay, this scenario also implies low D and 3 He and 
high Li, even if for a completely different reason (low iVjij? instead of high ?7bbn)- The 
neutron concentration is however affected in two opposite ways: higher initial n abundance 
and post-decay low iV e g giving neutrons more time to decay. These effects tend to compensate 
each other. The time that neutrons have to decay depends of how close to BBN the decays 
happen and indeed we find that 4 He grows as the ALP decay happens closer to BBN and 
a small 4 He region at ~ MeV and r ~ 30 s, where neutron decay plays a role inducing 
low He. However, because of the two opposing effects, the 4 He abundance is not a sensitive 
indicator of ALPs in this region. 

In order to numerically evaluate the impact of decaying ALPs in the BBN predictions, 
we have used a BBN code that includes the modified cosmology driven by ALP decays 
computed with the tools of the previous section. We have written a simple BBN code 5 in 
Mathematica to compute the primordial abundances of D, 3 He, 4 He, 7 Li and 7 Be. We have 
used the minimal reaction network relevant for rj ~ ?7cmb an d N c g ~ 3 as detailed for instance 
in [24, 25, 27]. This allows us to easily compute the outcome of BBN when ALPs have a non 
trivial role during BBN itself, and gives the right trend when ALPs decay much later than 
BBN, enhancing enormously the value of t/bbn with respect to r/cMB- 

The predictions for the primordial mass fraction of 4 He, Y p = 4nn e /nB, and the 
deuterium-over-proton ratio D/H in the ALP-decay scenario are shown in Fig. 5. The iso- 
contours very much resemble those of iV e gf because the outcome of BBN is mostly sensitive 
to the value of t/bbn and therefore to a possible baryon dilution, which qualitatively follows 
the same logic as the neutrino dilution. ALPs with small mass and fast decay disappear from 
the bath in LTE with photons at temperatures ~ m^. Therefore BBN only depends on the 
ALP mass, and not on their lifetime. When the lifetime is longer, the isocontours are parallel 
to the lines of constant entropy production as in the iV e ff case. We discuss the constraints at 
the end of this section. 

5 When facing standard ALP-less cosmology, our results are in very good agreement with standard BBN 
calculations obtained with the KAWANO [26] or PArthENoPE [27] codes, given the theoretical and experimental 
uncertainties, which gives us confidence in our results. 
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Figure 5. Isocontours of the primordial abundance of deuterium normalized to protons (D/H) (left) 
and helium Y p (right) in the decaying ALP cosmology, as a function of the ALP mass and lifetime. 



Large masses 

Photons from the decay of ALPs heavier than twice the photoionization threshold of elements 
can strongly modify the predictions of BBN. At high densities and temperatures, high-energy 
photons interact very fast with the thermal bath creating electromagnetic cascades, whose 
spectrum features a very sharp cut-off at high energies E > Eq ~ rrig/(22T) [28]. At high 
T, this cut off will lay below the photo-dissociation threshold of nuclei and these effects are 
negligible. We have followed the methods of [9] to check in which region of parameter space 
this phenomena can change our predictions. The photo-dissociation of deuterium (Q = 2.23 
MeV) will of course reduce its abundance further, making our bounds even stronger. On 
the other hand, the photo-dissociation of 4 He can create deuterium, which reverses the trend 
of our constraints. This requires higher-energy photons since the Q-value is 19.81 MeV and 
therefore a lower decay temperature. We have depicted the areas where these effects are 
relevant as grey regions in Fig. 5. 

If ?7j0 is larger than twice the charged pion mass m n + = 139.57 MeV, the decay channel 
4> —7- 77r + 7r~ opens up. Even if its branching ratio is very small, the abundance of relic ALPs 
is almost thermal and thus a huge amount of pions (compared with that of the present nuclei) 
can be produced. If the decay happens before or during BBN, the universe is dense enough so 
that pions can induce neutron-proton interconversions 7r + +n — > 7r°+p + and 7r~ ' +p + — > 7r°+n 
before decaying [30]. The second reaction is favoured because of the Sommerfeld enhancement 
and the typical overabundance of protons over neutrons at T < Q pn . These reactions will 
therefore tend to increase the neutron-to-proton ratio (they can do it much more drastically 
than the mere presence of the ALP during the p n freeze out commented in the previous 
subsection). The higher neutron abundance would aid the heavy element production but it 
also increases D/H. This is the most important effect. Since almost all neutrons end up in 
4 He nuclei taking protons with them, a higher initial neutron abundance yields a smaller 
final proton abundance and thus a larger D/H ratio. If n n /n p ~ 1 at the onset of BBN, all 
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protons end up in 4 He and D/H would be arbitrarily large! 6 . 

We have included the effects of pions in our BBN code following Ref. [31]. In Fig. 5 
we can see that the low D/H trend of low mass ALPs is drastically changed when the mass 
gets above 2m n + and the effect on 4 He gets strongly boosted when crossing this boundary. 
The effects of pions are hampered if the ALPs decay very early (r < 10~ 2 s), when the 
electroweak reactions p + + e~ -H- n + v e can still re-establish the n n /n p equilibrium; or very 
late (above r ~ 10 2 s), because pions fail to interact before decaying 7 . 

For ALP masses above a few GeV, ALP decays will produce quark-antiquark pairs that 
will hadronize. Hadronic cascades can dissociate nuclei, if they are happening after BBN 
(typical time of r ~ 10 2 s) [29, 31, 32]. Because of the large ALP relic abundance, the effects 
of electromagnetic or hadronic cascades are necessarily dramatic. We consider extremely 
unlikely that the combined effect of nonstandard BBN with the post-BBN processing gives 
similar results to standard BBN. So we exclude all the regions where electromagnetic and 
hadronic cascades play a role, see Fig. 5. In the literature, cascade constraints are usually 
presented in terms of m^n^/n^, plotted in function of lifetime r. In the case in exam, it is 
possible to refer to this representation almost directly from Fig. 5, considering that the ratio 
n^/n-y, given by Eq. 2.5, is constant in all the plotted parameter space in which cascades 
play a role and does not change much outside. 

Constraints 

We can now establish constraints by comparing our predictions with the measurements of 
primordial abundances present in the literature. 4 He was in the past the favourite indicator of 
the presence of extra degrees of freedom in the very early stages of BBN. The value of Y p can 
be estimated from an extrapolation to zero metallicity of the measured 4 He content of metal- 
poor extragalactic HII regions. The systematics of the measurements, the extrapolation to 
zero metallicity together with a somewhat unknown early stellar nucleosynthesis have caused 
the best estimate of Y p to vary significantly over the years, a fact that calls for extreme 
caution when quoting bounds. Here we shall be conservative and adhere to the proposal 
made in [36], where the authors set a robust upper bound on Y p based on the assumption 
that the helium content is an increasing function of the metallicity of the cloud. They find 

Y p < 0.2631 (95%C.L.). (3.12) 

The corresponding exclusion is depicted as a purple region in Fig. 6, the summary plot of 
this section, and in Figs. 1 and 2. 

When it comes to deuterium, we face other problems. The data are scarce, we only 
have reliable estimations from 7 high redshift low metallicity clouds absorbing the light of 
background quasars. The results of these estimations agree well at first glance but there is 
a scatter of the measurements beyond the expectations from the quoted systematics. The 
PDG quotes D/H= (2.82 ± 0.21) x 10~ 5 for the primordial ratio of deuterium to hydrogen, 
where the error has been enlarged to account for the still unexplained scatter. It is possible 
to obtain robust bounds from this measurement since the abundance of deuterium is known 
not to increase by the effect of unknown stellar processing. Any measurement of deuterium 
constitutes therefore a lower limit to the primordial value. In order to be conservative in this 

6 Of course including neutron decay during BBN would still give a finite, albeit very large, result. 
7 In this last tiny region our results cannot be taken quantitatively on trust, since we have not taken into 
account the possibly ineffective slowing down of pions after e + e~ annihilation. 
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Figure 6. Bounds on early ALP decays from deuterium underproduction (D, red), helium over- 
production (He, purple), helium photodissociation (7He, pink), hadronic cascades (Hadr, pink) and 
neutrino dilution (N e s, yellow). 

study we use 

D/H| p > 2.1 x 10" 5 . (3.13) 

The corresponding exclusion bound is shown in red in Figs. 6, 1 and 2. 

Note that the PDG average agrees nicely with the outcome of BBN calculations with 
the value of ?7cmb measured by WMAP, (we obtain D/H| p = (2.4 ± 0.1) x 10 5 where the 
error comes from the uncertainty in t/cmb) providing one of the most beautiful tests of 
standard cosmology. A word of caution is however in order since the WMAP value depends 
on cosmological priors such as the spectral index of primordial fluctuations. The quoted 
value stems on a scale-free power-law which we carry as a further assumption. 

The 7 Li abundance cannot be reliably used to constrain ALP decays in the region of 
interest since at the moment observations do not agree with standard BBN predictions. The 
general trend is to increase the discrepancy between predictions and observations, since 7 Li/H 
increases with ??bbNi with a low iV e gf and a high n n /n p at BBN. This observation disfavours 
any reasonable attempt of solving the Li problem with decaying ALPs. Awaiting the Li 
problem to be solved, we have chosen not to propose any additional constraint although 
we observe that might theory and observations be reconciled, the sensitivity would be very 
similar to that of deuterium. 

3.3 Summary 

Our summary plot for bounds on early ALP decays is Fig. 6. Deuterium, helium or N e g 
can be the most constraining argument depending on the ALP mass. Let us comment on 
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why this is so. For very small lifetimes, ALPs will be in thermal equilibrium and the bounds 
only depend on the ALP mass. In this case deuterium dominates. The bound corresponds 
to ALPs that inject a a fraction of order 10% of the total entropy in the electromagnetic 
bath after BBN. Under this circumstances the baryon and neutrino dilution is 0(1), and, 
since deuterium is the most sensitive observable to t?bbn> it is also the most constraining 
argument. Note that the 4 He abundance depends only logarithmically on tjbbn while D/H 
oc 1/(?7bbn)~ 1 ' 6 - Moreover, we have already said that the N e g bound is not enough to 
constrain the LTE decay of ALPs. 

At larger masses, ALPs in LTE would disappear before BBN and we loose the constraints 
for arbitrarily small lifetime. However, as the ALP lifetime increases, the decay proceeds 
increasingly out of equilibrium together with the entropy injection and the neutrino and 
baryon dilutions. For uia, < 10 MeV we can only constrain decays happening during/after 
BBN (r > 10 2 s) for which D/H is still the most constraining argument. But as the ALP 
mass increases, so it does the dilution and for > 10 MeV it is big enough to be constrained 
by N e g. The bound on r improves quadratically with (for a fixed entropy injection we 
have m^Y^y/r = const., c.f. Eq. 3.9) until lifetimes of the order of the neutrino decoupling 
and then flattens. The deuterium bound follows quite closely behind N e g because a low N e g 
during BBN also implies a low D/H. Note that for lifetimes smaller than 10 2 s, i]bbn = ^cmb 
so the baryon dilution does not play any role and the low Neg is the only responsible for the 
low D/H. 

For ALP masses above 2m^+, the few pions produced in ALP decays affect notably 
the neutron/proton equilibrium before BBN leading to unacceptable high 4 He. This effect 
dominates over the neutrino dilution because a very small amount of pions is needed to change 
n n /n p while entropy injection must be 0(1) to significantly affect N e g (there are 10 9 times 
less nucleons than photons!). However, the N e g bound, even relying on very different physics, 
lies very close. This can be easily understood. In this region the entropy injection is huge and 
the ALP rest mass dominates the energy budget of the universe at the decay. Nevertheless, 
the effects we are constraining are 0(1) changes to standard cosmology. Therefore our 
constraints correspond to cases where most of the entropy release has been absorbed by 
thermalization processes, i.e. most of the decays happen when the universe is strongly secured 
against non-thermal distortions. However, the very last ALPs, whose presence does not even 
affect the expansion significantly, are still able to produce an observable effect. The number 
of these ALPs depends exponentially on r and therefore any bound on r can only depend 
logarithmically upon any of the other quantities of the problem, such as the ALP mass, 
hadronic branching ratio, etc. 

Finally note that all bounds disappear for fast-decaying large-mass ALPs. These ALPs 
disappear from cosmology before the freeze-out ofpffn reactions and neutrinos leaving no 
trace in the output of BBN. Here the ALP mass is crucial, because fast-decaying low mass 
{ m 4> ^MeV) ALPs are kept in equilibrium with photons through the inverse decay processes 
77 — > (f) down to temperatures where they affect BBN, even if their lifetime would have made 
them decay much earlier. This is the reason why the BBN bounds do not disappear at low 
masses even for very short lifetimes in the figures of this section. However, in this section 
we have assumed the decay of ALPs to occur before the CMB release epoch, in order to 
trust the standard analysis of CMB data that give us values for ijcmb and N e g. Therefore, 
we can not formally use our results to constrain ALPs which are kept in equilibrium longer 
than the time at which N e g and r\ are imprinted on the CMB. To be conservative, we take as 
a boundary time the matter-radiation equality, whose temperature is T^q ~ 3 gVj and. thus 
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our bounds cease at tua, ~ 3T cq ~ 10 eV, see Fig. 2. Anyway, the existence of ALPs which 
decay in LTE and have a mass smaller than 10 eV is unconvincing from the cosmology point 
of view, as it would affect the CMB, for instance altering the transparency of the universe 
to photons. Besides this, such high coupling and light mass ALPs are definitely excluded 
by stellar evolution, positronium and T decays and helioscope searches [3], which makes less 
important to know exactly where our bounds are no more valid. 



4 Late ALP decays 

We have just shown how it is possible to exclude some ALP parameter space if the decay 
happens before matter-radiation equivalence. Later decays can also be constrained, directly 
measuring the relic photons emitted or checking the distortion on the CMB spectrum they 
would create. 



4.1 Direct detection of ALP decays 

After recombination the universe becomes practically transparent to radiation, since almost 
all the electrons are captured by nuclei forming neutral atoms. The photons injected by ALP 
decay can be in principle directly detected, unless their wavelength lies in the ultraviolet 
range and they are absorbed in the photoionization process of atoms. 

In the parameter space we can constrain, the ALP decays happen at rest in the comoving 
frame. The spectral flux of photons produced in the decay of a diffuse ALP population 
is [11, 37] 
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where the subscript means quantities at present time, n^o is the putative ALP number 
density if <f> would be stable and Eq is the energy at which the photon would be seen today. 
The photon initial energy is m^/2 and thus the redshift of the decay is given by 1 + z = 
(t7u/2)/-Eo, i-e. the ratio between the emitted (then) and the measured (today) energy. For 
simplicity we assumed matter domination neglecting 0(1) corrections due to the cosmological 
constant. 

Radiation above 13.6 eV can photoionize hydrogen, an effect that we take into account 
correcting the flux by multiplying Eq. 4.1 with the survival probability 

P( z ) = e ~^ E ) (4.2) 
f z dz' 

<z,E) = J ^————n H (z')a pe (E), (4.3) 

where E = Eq(1 + z), nn is the hydrogen number density and a pe = 256n(Ei s /E)^ / (3am^,) 
is the Hydrogen photoelectric cross-section with E± s = 13.6 eV. We have thus compared 
the resulting spectrum with the extragalactic background light (EBL) spectrum we have 
found in the review by Overduin and Wesson [38], reproduced in Fig. 7. The region of 
parameter space excluded by this comparison is plotted in dull green in Figs. 1 and 2 and 
it is labelled EBL. The same approach was used in Refs. [38, 39] to look for axions in the 
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Figure 7. Extragalactic photon background [38]. Courtesy from J. Overduin and P. Wesson. 



optical EBL. In this parameter space the ALP abundance would be the minimum considered 
n^/riry = (3.36/106. 75)/2 ~ 0.016. As we noted, for > 154 eV the ALP energy density 
overdoses the universe, in gross contradiction with observations. Therefore, above this mass 
we present our bounds assuming that ALPs provide the right amount of DM, which of course 
requires some non-standard dilution of the ALP number density by additional degrees-of- 
freedom above the electroweak scale. 

The bound gets severely degraded in the = 13.6 — 300 eV range where not only 
absorption is very strong but also the experimental data are extremely challenging and the 
EBL spectrum has only an upper bound estimate (of course these facts are closely related). 

A more powerful strategy to search for decay photons is to examine the light emitted in 
galaxies and large scale structure, where the dark matter density is above the average. Also 
the backgrounds are above the EBL, but the decay photons show up as a peak at frequency 
u> = in the galactic spectra, which helps background subtraction. Searches in the 

visible have been presented in Refs. [39, 40] and more recently in Ref. [41], and in X-rays in 
Refs. [42, 43] searching for sterile neutrino dark matter decays. We followed these references, 
rescaling their results for the ALP decay mode and abundance (for < 154 eV), and we 
obtained the exclusion bounds in Figs. 1 and 2 labelled respectively Optical and X-Rays. 
We have also plotted in Fig. 2 the exclusion bound obtained by galactic line searches in the 
7-ray range for decaying DM, labelling it 7-rays. We took the data from [44] and again we 
rescaled them for the ALP decay mode. 

The radiation produced by ALP decay, being proportional to the relic ALP density, 
would be linearly sensitive to the particle content above the EW scale through <?*(Tf ). The 
constraint on g, for < 154 eV, relaxes as yl + new dof/106.75 with the new thermal 
degrees of freedom with masses between the EW and Tf Q . 
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4.2 Distortions of the CMB spectrum 

The constraint that allows to fill the gap between the BBN and iV e g excluded region and 
the telescope bounds (and also largely overlaps with these) is given by distortions of the 
CMB spectrum. FIRAS measured the CMB spectrum to be a perfect black-body within the 
experimental errors [45, 46]. Any amount of electromagnetic radiation from the ALP decays 
should have had enough time to fully thermalize or must be very tiny compared with the 
CMB itself, O(10~ 5 ) [45]. The ALP energy injection is always greater than 10~ 5 , and thus is 
potentially in radical conflict with observations, unless ALPs have not decayed significantly 
yet. On the other hand if the decay is very early, the distortions can be r ether malized. 
Evading these arguments requires cosmological scenarios with g*(Tfo) > 0(1O 5 ). 

The most important processes for the thermalization of injected photons are double 
Compton scattering, e77, and Bremsstrahlung pe <-> pej, which permit to fully re- 

gain the thermodynamical equilibrium distribution. In standard cosmology double Compton 
dominates 8 because of the small value of i], being extremely effective until the temperature 
drops below Tdc ~ 750 eV. This leaves Compton scattering cy — > as the most relevant 
interaction between photons and the rest of the plasma. Through Compton scattering only 
kinetic equilibrium can be obtained, as it can redistribute the energy of photons but not 
change the photon number, thus leaving a non-zero degeneracy parameter [47] 



J_ (qP4> 
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(4.4) 



T=T d 

which is constrained by FIRAS to be |^| < 0.9 x 10 [45]. ALPs with decoupling temperature 
smaller than Tdc are excluded while in the opposite case the produced distortions can be 
efficiently erased. The light green region labelled CMB /i in Figs. 1 and 2 corresponds to 
Td < Tdc- Note that this bound is somewhat conservative, as a significant amount of energy 
can be released after T<j, especially above ~ keV where the ALP energy dominates the 
universe before the decay. 

At Tq ~ 25 eV also Compton scattering freezes-out and the above bound has to be 
reformulated. Electrons rapidly thermalize with the non-thermal population of photons, but 
the CMB can not drain efficiently energy out of them. The degree of thermalization that 
CMB photons can still attain depends on their energy and this imprints a typical pattern 
on the CMB spectrum. The parameter y, characterizing this distortion, can be estimated 
as [11] 



Pi 



~ exp(4y) - 1 (4.5) 

T=T d 



and is bound to be |y| < 1.5 x 10 5 [45]. The region excluded by this last bound is roughly 
Td < Tq and is shown in Figs. 1 and 2 coloured in light green and labelled CMB y. 

4.3 Reionization history 

At temperatures around 0.3 eV (redshift z ~ 1300), most electrons and protons combine into 
neutral Hydrogen and the universe becomes effectively transparent to the CMB. The universe 



Actually in the decaying ALP cosmology this will be the case as well. Even if the value of r\ before the 
decay can be much lager than -qcMB, the amount of photons injected by ALPs is so huge that even the latest 
ALP decays are dangerous. The latest ALPs are not energetically significant and decay when r\ ~ »?cmb so 
double Compton will dominate in their thermalization. 
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reionizes again much later, between redshifts 6 ~ 10 presumably due to ultraviolet emission 
from the first galaxies but the details of this process are still poorly understood. The details 
of the reionization process can be studied by the slight imprint that the free electrons leave 
in the CMB through Thomson scattering, for instance in the polarization. The optical depth 
for CMB photons is one of the parameter that can be measured from the CMB and it is 
defined as 

r opt (z 1 ,z 2 ) = -j^ H{z){l + Z) dz, (4.6) 

where cry is the Thomson cross section, n e is the total electron density and Xi on the fraction of 
them not trapped in nuclei, i.e. the free electron or ionized fraction. The WMAP7 measured 
T op t after recombination to be 0.088 ±0.015 [23]. A factor 0.04 — 0.05 of this can be attributed 
to a fully ionized universe up to redshift ~ 6, as supported by the absence of Ly-a features 
in quasar spectra. The origin of the remaining fraction, rg ~ 0.04, is still uncertain. 

We already commented that it seems hopeless to detect directly ALPs decaying into 
ultraviolet photons (energy range 13.6 ~ 300 eV), because they photoionize very efficiently 
Hydrogen and they are thus absorbed very fast. Of course, these ALPs contribute to the 
reionization story and they can be constrained indirectly, through their contribution to the 
CMB optical depth. 

It is simple to compute a rough estimate of the ionization fraction x; on caused by the 
ALP decay photons [48] . Assuming that each decay photon ionizes only one H (immediately 
after the ALP decay) the number of ionizations per unit time can be estimated as 

~ 2r 7 ^- ~ 2 x HT 3 (-^-Y ( — ^ Xe-'A Myr" 1 . (4.7) 

^ W 7 n^(z) VlOOeW VlO" 13 GeV~V 

If we now multiply this quantity with a typical time scale [48] 

t H = l/H(z) ~ 2.4 Myr (501/(1 + z) f /2 , (4.8) 

we get a conservative estimate of the typical degree of ionization induced by the ALP decays 
up to a certain redshift. ALPs with 100 eV mass and g ~ 10 _13 GeV _1 would have produced 
an ionization comparable with the standard residual value 4 ~ 10~ 4 already at high redshifts 
~ 500, and the ionization grows in time as (l + z) _3//2 showing a potentially interesting effect. 
Indeed one could even think that ALPs close to these parameters may have had a role in 
reionizing the universe, but the (1 + z) -3 / 2 dependence is too soft. Reionization seems to be 
a much more abrupt process. 

In order to obtain a more detailed constrain, we calculated the ionization history of 
the universe in the decaying ALP cosmology by introducing the ALP ionizations in the 
recombination code RECFAST [49]. Then, we computed the optical depth in the interval z = 
6—100, requiring it to not exceed tq. In the calculation, we used the ALP thermal abundance 
and considered only a ionization for each decaying ALP, which is certainly conservative. Our 
results are excluding the light green region labelled Xi on in Figs. 1 and 2. This bound would 
increase up to one order of magnitude at the largest masses for which ionization is effective, 
m 4> ^5 300 eV, if we assume that all the energy of the emitted photons can be converted into 
ionization. 
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5 Beyond the two photon coupling 



A generic pseudo Nambu-Goldstone boson can feature other couplings to SM particles [50]: 
anomalous couplings to two gluons 

£ 3 c 39 ^^Tr{G^G^} (5.1) 

and derivative couplings to fermions /, possibly even flavour non-diagonal 

C3c ff Mf l5 ^f (5.2) 

with, in principle, 0(1) coefficients c gg , Cff. The addition of other couplings to SM particles 
necessarily implies a smaller temperature for the ALP decoupling and therefore a higher 
primordial abundance. In this sense, considering only the two-photon coupling for the relic 
ALP production produces conservative bounds. All what remains is the ALP decay, for 
which the new couplings above open new decay channels. The coupling to fermions allows 
the decay 4> — > ff at a rate 

\ 2 / , \ 2 / / T~2 / / \ 2\ 3 / 2 

Cff \ (mf + m,f>) / fmf + mf\ I jmf — mp 




r ^-[t) ar Y -r^rJ I'-l^ 1 1 • (53) 

which is suppressed with respect to the two-photon decay for small fermion masses. Writing 
g = c 77 a/(27r/^) the (f> — > ff can dominate only in an interval near the kinematic threshold 
1 > (rrif + m'f)/m0 > ac-y-y/Ancf f . For ALP masses above few GeV, the coupling to gluons 
allows the ALP decay into two gluons at a rate T^ gg = 8(c 33 /c 77 ) 2 r^^. 77 . At low masses 
this coupling implies a phenomenology very similar to the axion case 9 . The ALP will mix with 
the rf and, through it, with the pseudoscalar mesons and hadrons and get new contributions 
to the two-photon coupling. 

Let us now review the impact of these new decay channels in our bounds. First of 
all note that if the ALPs are cosmologically stable, the bounds from direct detection of 
ALP decay photons and the DM overproduction still hold. These arguments span the long 
lifetime range of our constrained region. The short lifetime region corresponds to BBN and 
N c q constrains and it is summarised in Fig. 6 in the ALP mass and lifetime plane. At low 
masses, the deuterium and He bounds come from ALPs in thermal equilibrium with the bath. 
Clearly, adding more couplings between the ALP and SM particles we cannot avoid these 
bounds. In the intermediate mass region 300 keV < < 2 ffin, where = 105.7 MeV is 
the muon mass, the D/H and N e Q bounds follow from the dilution of baryons and neutrinos 
with respect to photons. These bounds apply to ALPs decaying into photons or electrons 
(actually we have not made a difference between the two in our equations) and since the 
direct decay into neutrinos is suppressed by ~ (ra^/m^) 2 , an amazingly tiny number 10 , the 
bounds are perfectly valid provided one interprets r as the total lifetime (not only due to the 



9 This particle cannot solve the strong- CP problem unless it is built-in massless or nearly massless. All our 
bounds apply to particles more massive than the axion so there is no point in considering this possibility here. 

10 Unless one considers sterile neutrinos with m„ ~ but then neutrinos have a strong tendency to 
constitute too much DM. A way to avoid this is to make them decay into SM neutrino + photon, but this 
produces entropy so we expect a similar, slightly smaller, bound from D/H in this case. In these models the 
low N e s tendency is reversed since the sterile neutrinos produce neutrinos in its decay. 
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two photon decay channel). If the decay into two electrons dominates, when we translate the 
bounds of Fig. 6 in the — g plane they will show worst than if we only consider the two 
photon coupling. The lower bound on g reduces by a factor ~ 2Trc ee m e /ac Ty m ( j ) . 

If 2m w + > > 2m^ we have a somehow different scenario where the ALP tends to 
favour the <f> — > /J + /^~ decay. The upper limit on r in this region comes from having too low 
N e g already before BBN. But If the decay into muons dominates we will rather have a high 
N e Q because the amounts of energy released in electrons and in neutrinos by muon decay 
[i — > evv are similar. Since data favours values larger than the standard N e g = 3 the N e g 
bound will relax somehow. We do not expect them it to disappear because ALPs can still 
produce too many neutrinos. Also in this case the bound on deuterium should come from a 
too high D/H, which is less conservative a constraint. In any case the bound from He will 
stay since it mainly comes from a high 7/bbn and the ALP contribution to the expansion at 
the freeze out of p n weak reactions. 

Finally, for > 2m n + the most stringent bound comes from 4 He overproduction due 
to the presence of charged pions before BBN enhancing the neutron/proton ratio. As we 
commented this bound does depend very little of the details and branching ratios of the ALP 
since only a minimal number of pions would do the job. Therefore we expect it not to change 
very much. However, then quoting this constraint in the — g plane this bound would 
display lower than in the case where only the two photon case is considered. Only in this 
region the coupling to two gluons can affect the ALP decay and will certainly increase the 
pion multiplicity of the decay making the bound on r slightly better. The decay into muons 
can dominate if is not too far from 2m ^ and all said in the above paragraph holds. It 
appears that the helium bound will still be the most relevant in this case. 

6 Conclusions 

In this paper we have evaluated the impact of pseudo Nambu Goldstone bosons featuring 
a coupling to two photons, usually called axion-like particles (ALP) and in particular those 
more massive than the QCD axion. This paper reviews and complements the previous work 
of Masso and Toldra [10, 11]. 

ALPs are efficiently created in the early universe via the Primakoff effect and then can 
decay into two photons leaving traces in the density of neutrinos, the primordial abundances 
of light nuclei (BBN), the spectrum of the cosmic microwave background or simply creating 
a diffuse photon background. We have found that the more stringent constraints for early 
decays (before recombination) are set by the density of neutrinos (A^g) and the primordial 
abundances of deuterium and helium, depending on the ALP mass. These bounds are sum- 
marised in Fig. 6. The diffuse photon background and the CMB spectrum are the most 
relevant for longer lifetimes. Interestingly, we have found that these bounds are only slightly 
modified in scenarios where the ALP has other couplings and other possible decay channels. 

At small masses and large couplings, these bounds are complemented by stellar evo- 
lution arguments, particularly the ratio of red giants to horizontal branch stars in globular 
clusters, which we have revisited to define precisely the high mass frontier. Altogether, these 
arguments exclude a huge patch of parameter space, shown in our summary plots Figs. 1 
and 2. 
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A The globular cluster bound for high mass ALPs 

The ratio of red giants (RG) to horizontal branch (HB) stars in globular clusters is an 
indicator of the presence of weakly coupled low mass particles which can efficiently drain 
energy from the stellar cores [51]. The non standard energy loss prolongs the RG phase and 
shortens the HB in a different way, being these effects sensitive to the type of particle and 
coupling to matter. For particles coupled to two photons, the lifetime of HB stars is more 
affected than the RG. Qualitative studies [51, 54] show that the exotic energy emitted per 
unit time and mass, averaged over a typical HB core, has to satisfy 

(e CX otic) < 10 g- 1 ergs" 1 (A.l) 

The main reaction producing ALPs in a HB core through the two photon coupling is 
the Primakoff conversion in the Coulomb field of ions in the plasma, i.e. 75 — > <pq where the 
qs are mostly He 2+ and protons. The cross section in the ion rest frame can be written as 

m 2 \ (2uj 2 (1 + f3) + k 2 -m 2 ' 




2uj 2 J b \2u 2 {l- ^) + k 2 
m*+k 2 (2u 2 (l + (3)-m 2 




(A.3) 



where j3 = y uj 2 — m^/u is the ALP velocity, we have assumed the ALP mass much smaller 
than the mass of the ion, so that the ALP energy is equal to the initial photon energy uj, and 
we have taken into account Coulomb screening following [55]. In an non-degenerate plasma, 
the screening scale is given by the Debye-Hueckel formula 

72 47ra ,. „. 

k 2 = —n q (A.4) 

where the sum involves all charged species in the plasma. At low masses, the logarithmic 
divergence of the cross section (inherited from the well known forward divergence of Coulomb 
scattering) is cut-off by screening. However, when considering large masses, it is the ALP 
mass which plays this role and screening becomes irrelevant. 

The energy loss per unit mass is given by folding the cross section with the photon 
phase space and thermal phase space distribution /be(^) = (e w / T - l)" 1 

e = ~ J (^3 w /be(w) Qjiija^uj) . (A.5) 

For typical conditions in a HB core, T ~ 8.6 keV and matter density p = 10 4 g cm -3 , the 
screening scale is k s ~ 27 keV. 
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Figure 8. Bounds on ALPs parameter space from energy loss in Horizontal Branch stars. 



The energy loss grows with g and so does the impact on the star. However, at some 
point g can be so large that the inverse Primakoff process, 4>q — > 75, cannot be neglected 
and ALPs will be reabsorbed inside the star. This does not invalidate the constrain, since 
ALPs will nevertheless transfer energy from the core to the external shells of the star which 
has a very similar effect on the star evolution [53]. However, if the ALP mean free path, A, 
gets very small, even the energy transfer will become eventually negligible. Unfortunately, 
the energy transfer on a HB core is convective and therefore it is difficult to ascertain for 
which exact value of g ALPs are harmless. In order to get an estimate, the authors of [53] 
propose to compare the ALP energy transfer with the radiative energy transfer mediated 
by photons. Since radiative transfer is subleading in HB cores, the limit of validity of the 
bound is conservative. Thus, we impose that the contribution of ALPs to the Rooseland 
mean opacity is smaller than the standard value in a HB core kq ~ 0.5 cm 2 /g, which means 

where / BE = <9/be/<9T, the mean free path is given by A" 1 = Zjrija^zi^) with the 
inverse Primakoff cross section given by detail balance a^zi^) — °7z(w)/(2/3 2 ). The full 
range excluded by HB stars in globular clusters is shown in Fig. 8 and is also reproduced 
Figs. 1 and 2. 
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